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We demonstrate control of the thermal hysteresis in superconducting constrictions by adding a 
resistive shunt. In order to prevent thermal relaxation oscillations, the shunt resistor is placed in close 
vicinity of the constriction, making the inductive current-switching time smaller than the thermal 
equilibration time. We investigate the current-volt age characteristics of the same constriction with 
and without the shunt-resistor. The widening of the hysteresis-free temperature range is explained 
on the basis of a simple model. 


INTRODUCTION 

A superconducting weak-link (WL), such as a constric¬ 
tion, between two bulk superconductors is of interest for 
its Josephson junction-like properties and subsequent ap¬ 
plication to micron size superconducting quantum inter¬ 
ference devices (/r-SQUIDs) [U [2]. The latter can be 
used in probing magnetism at small scales [3H6] . Hys¬ 
teresis present in current-voltage characteristics (IVCs) 
is a limiting factor in WL-based SQUIDs. In a hysteretic 
IVC when the current is ramped up from zero, the device 
typically switches to a non-zero voltage state at the crit¬ 
ical current /c- The subsequent current ramp-down gives 
a switching to zero-voltage state at a smaller current, 
called re-trapping current /^. Hysteresis in IVCs is seen 
at low temperatures and disappears above a crossover 
temperature Ty as Ic and Ir meet [zm]. In a conventional 
tunnel-barrier type Josephson-junction, hysteresis arises 
from large junction capacitance and can be eliminated by 
adding a shunt-resistor in parallel to the junction mm- 
The effect of the shunt resistor on nano-wire based WL 
devices was modeled recently using resistively and capac- 
itively shunted junction (RCSJ) model with an effective 
capacitive time m- The hysteresis in similar devices 
is well understood using the thermal model [12]. The 
hysteresis in WLs is due to local Joule-heating [13114], 
which gives rise to a self-sustained resistive hot-spot in 
the WL region, even below Jc- 

Eliminating thermal hysteresis in WLs has been the 
subject of intense research in the past years. A normal 
metal shunt directly on top of the constriction [HHni has 
been tried, but it affects both the superconductivity and 
thermal properties in a way that depends on the inter¬ 
face transparency. Using a bilayer with a superconductor 
(that can locally etched with a Focussed Ion Beam) cov¬ 
ering a normal metal film allows one to obtain a WL that 
is also a good thermal bath m- A parallel shunt resis¬ 
tor far away from the WL m is more flexible approach, 
but it gives rise to relaxation oscillations due to large 
inductive time for switching of the current between the 
WL and the shunt. The performance of such SQUIDs 



FIG. 1: (a) SEM image of the WL device after etching the 
gold shunt. The inset in (a) shows the zoomed-in SEM image 
of the WL with patterned width 70 nm and length 150 nm. (b) 
SEM image of the same device rotated by 90°, before etching 
the hatch pattern shunt resistance made of gold microwires, 
(c) Resistance vs temperature (at 0.01 mA bias current) for 
the shunted and unshunted WLs, showing different transitions 
at Tci and Tc 2 respectively corresponding to the narrow and 
the wide leads. 


with a distant shunt-resistor is eventually similar to that 
of the hysteretic ones mm- A systematic study of the 
ability of a parallel shunt in preventing both the thermal 
runaway and hysteresis is thus highly desirable. 

The role of a shunt-resistor can be understood using a 
simple quasi-static thermal model discussed by Tinkham 
et al. m- In this model, the heat generated in the re¬ 
sistive hot-spot in a long constriction is conducted (only) 
through electronic conduction to the large electrodes at 
the end. The thermal conductivity K of the normal metal 
and superconductor are assumed to be identical and 
temperature-independent. The re-trapping current is 
then found to be Ir{Tb) = Ir — T})/Tc with /r(0) = 

^^KATr/LRn. Here Rn is the normal resistance of the 
constriction of length L and cross-sectional area A, is 
the bath temperature and is the superconductor crit- 
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ical temperature. From Ginzburg-Landau theory [2] the 
critical current follows IdT^) = /c(0)(l — T^h/Tc) in the 
regime > Tc/2. Thus /c and cross at a crossover 
temperature = Tc[l — [/^(0)//c(0)]^] [20]. In the pres¬ 
ence of a shunt resistor Rg, the bias current is shared 
between the shunt and the WL, when the latter is re¬ 
sistive. Thus, in /^(O) expression, 1/Rn is replaced by 
(l/i?n) + (1/^s)- As a result /r(0) changes to a higher 
value given by 


/^,(0) = Iri0)^/1 + Rn/Rs. (1) 

In contrast, Ic remains unaffected. Hence the crossover 
temperature decreases and the hysteresis-free tempera¬ 
ture range [Th^Tc] widens thanks to the shunt. For elim¬ 
inating hysteresis above temperature T, Rg with value 
less than Rsc = Rn/[{Ic{T) / Ir{T)d — 1] will be required. 
The assumed immediate sharing of the bias current be¬ 
tween WL and Rg implies a small inductive current¬ 
switching time as compared to the thermal equilibration 
time. The minimum shunt resistor value, Rgc^ from our 
simple model can also describe the behavior of shunted 
nano-wire devices studied by Brenner et al. m- 

In this Letter, we compare the current-voltage char¬ 
acteristics of carefully designed WL devices with (and 
without) a shunt resistor kept in close vicinity of the WL, 
thus making the inductive current-switching time smaller 
than the thermal equilibration time. We observe an in¬ 
crease in the re-trapping current and a widening of the 
hysteresis-free temperature range thanks to the shunt, 
which we discuss using the above model. 


EXPERIMENTAL DETAILS 

Devices were fabricated on Si substrates in two sub¬ 
sequent lithography and e-beam deposition steps as fol¬ 
lows: 1) laser lithography of hatch patterned shunt re¬ 
sistor [21] and alignment marks on a photo-resist, 2) de¬ 
position and lift-off of a Ti (3 nm)/Au (20 nm) layer, 
3) oxygen reactive-ion-etching (RIE) to remove residual 
resist, 4) deposition of a 31 nm thick Nb-film, 5) aligned 
electron beam lithography of a PMMA resist of the WL 
pattern, 6) deposition and lift-off of a 20 nm thick Al- 
film, 7) etch of Nb with SFe-RIE, 8) chemical removal of 
Al. After fully characterizing the shunted device, the An 
shunt was etched using a KI-I 2 solution which does not 
attack Niobium. Electrical transport studies down to 1.3 
K were pursued using a closed cycle He-refrigerator [22] 
with a homemade sample holder that incorporates cop¬ 
per powder filters. The data were recorded using data 
acquisition cards and homemade analog electronics. Two 
nominally identical devices demonstrated similar results. 



FIG. 2: (a) I-V characteristics for the unshunted weak link at 
a set of different temperatures. 


RESULTS AND DISCUSSIONS 

Fig. [^ shows the electron micrograph and resistance vs 
temperature of the reported device. The WL [see inset 
of Fig. [^a)] as designed has a 150 nm length and a 
70 nm width. Narrow leads (width 0.3 /im, length 2.4 
jiim) with normal resistance 2Ri, connect symmetrically 
to the two ends of the WL of normal resistance Rwl- 
Wide leads of width 2 fim connect the narrow leads to 
the shunt and the voltage probes. Fig. [^c) shows the 
measured resistance R versus temperature T before [see 
Fig. [^b)] and after [see Fig. [^a)] etching the An shunt. 
For the two cases, R drops from a saturation value of 
135 or 72 U (at 10 K) with two transitions at Td = 

8.5 K and Tc 2 = 8.8 K. These two critical temperatures 
are attributed respectively to the narrow-leads and the 
wide leads. Let us stress that the critical temperatures 
Tci ,2 are not affected by the removal of the An shunt, 
which confirms that the shunt etching did not damage 
the niobium pattern. From the resistance drop at Tc 2 for 
the unshunted device, we find a square resistance R^ = 

3.5 U giving a resistivity value of 10.8 fiQ.cm for the Nb 
film. From the narrow lead resistance and dimensions, we 
estimate Rwl = 9 U and 2Ri = 56 U. By comparing the 
resistances of the two devices just below Tc 2 , we find Rg = 
35.6 U, which is consistent with separate measurements 
of An wires’ resistances. 

While designing the device, we have kept the shunt- 
resistor close to the WL to minimize the associated loop- 
inductance Lgh. The inductance of a square loop with 
sides a and width 6 ^ a, is given by ^/j.oa\n[^] [23] . 
Using this relation with a = 20 jam and b = 3 /im, we 
estimate the inductance of the loop containing gold shunt 
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FIG. 3: (a) I-V characteristics for the same weak link device 
with a parallel shunt, at the same set of different temperatures 
as in Fig. 


hatch pattern and the Nb leads as Lsh = 40 pH. This 
gives a inductive-current-switching-time tl = Lsh/{Lis T 
Rwl) — 1 ps. The heat is transferred to the substrate 
over a length-scale given by thermal healing length lth{= 
^JKt/a = 1.6/im) j8| [24]. Here t is the film thickness 
and a is the interface heat loss coefficient. Thus the 
thermal cooling time is tt = Ith^/' k^D = ct/'K^a with 
D {=K/c = 1 cm^/s) as the diffusion constant and c 
as the volumetric heat capacity. We thus estimate the 
thermal time tt = 2.5 ns, which is much larger than 
Tl. When the WL switches from the superconducting to 
the resistive state, the current redistribution between the 
shunt and the WL is thus much faster than the thermal 
runaway in the device. 

Figureshows IVCs of the unshunted device at various 
temperatures. At low temperatures, sharp voltage jumps 
and drops are observed at the critical current Ic and at 
two re-trapping currents Ri and Ir 2 - The latter two arise 
from thermal instabilities respectively in the WL plus 
the narrow leads (/ri) and in the wide leads (/r 2 ) [H]- 
At 1.3 K, Ic is higher than both Ri and A 2 [see Fig. 
[^a)]. The IVC slope above 1^2 is 142 11, slightly larger 
than the normal state value of 135 Q because of over¬ 
heating. At higher temperatures when A 2 is smaller and 
heating less, the slope is 135 H. The IVC slope above 
Iri is 73 H, which is close to the combined resistance 
Rwl + 2i?i, i.e. 65 H. The slightly larger value is due 
to the spread of the hot-spot into the wide leads. With 
increasing temperature, R crosses A 2 near 3.2 K [see Fig. 
2 [b)] and it merges with Ri near = 6.25 K [see Fig. 
2[e)]. At higher temperature, the IVC is non-hysteretic 
and the resistance for /> Ic is 65 H, indicating that the 


WL as well as the narrow leads are resistive. 

Figure shows IVCs of the same device but prior to 
the shunt removal. We observe voltage jumps and drops 
at R and Ri while the second retrapping current R 2 
is visible only in IVCs with a larger bias current excur¬ 
sion m- In the resistive region, the slope is always 22 
Q which corresponds to the parallel combination of the 
normal resistance of the WL plus the narrow-leads with 
the shunt, i.e. {Rj^ + (2i?i + Rwl)~^)~^ • The critical 
current R magnitude at low temperatures is the same, 
within the error bars, as that of the unshunted device, 
confirming that the shunt removal did not damage the 
WL. Remarkably, the re-trapping current Ri has a higher 
value as compared to that of the unshunted device. As 
a result of Ri enhancement, R and Ri meet at a lower 
crossover temperature T^s = 5 K in the shunted device, 
see Fig. |^d). 

We summarize the temperature dependence of R^ Ri 
and R 2 for both devices in Fig. In every case, the 
retrapping current R 2 nearly follows a square root de¬ 
pendence with the bath temperature [8] extrapolating to 
zero at Tc 2 . This is consistent with R 2 being related to 
the thermal instability of the wide leads. Between about 
3 K and the crossover temperature or T^s , the critical 
current R of both devices decreases linearly with tem¬ 
perature. The extrapolated critical temperature R close 
to 7.2 K is that of the WL itself. In both the devices, 
the critical current R decays markedly slower above R 
or Rs , owing to the proximity effect [8] . 

The retrapping current Ri follows a similar temper¬ 
ature dependence in the two devices but with a higher 
magnitude, by a factor of about 1.64, in the shunted de¬ 
vice. This factor is similar to that found from Eq. 
i.e. 1.67, by using Rs = 35.6 O and Rn = 2Ri -f Rwl = 
65 11. A more appropriate model, incorporating the in¬ 
terface heat loss, for our device configuration provides a 
similar agreement, where the ratio is found to be close to 
1.60 [24]. For our devices, both the models give similar 
agreement 1211, as the conduction dominates over inter¬ 
face heat loss because the narrow leads’ length is compa¬ 
rable to Ith |H| • The two models have significant disagree¬ 
ment when the constriction is much longer than R [24|. 
Thanks to the shunt and the related Ri enhancement, 
the hysteresis-free temperature range has increased from 
[6.25 K, 8.6 K] to [5 K, 8.6 K], see Fig. For instance, 
the WL without shunt is hysteretic at 5 K, see Fig. 2 [d), 
while the one with shunt is non-hysteretic, see Fig. 3 [d). 

The merging of R and Ri above R is different from 
our earlier results on unshunted /i-SQUIDs |8| where a 
crossing of the two was seen at R. Due to the presence 
of the SQUID loop, the heat evacuation in the /i-SQUIDs 
is more efficient. For single WL devices with similar R 
values, the less efficient heat evacuation favors merging 
over crossing. In fact, just below R^ where R is small, 
heat evacuation eventually dominates, and we do see dis¬ 
tinct signatures of both R and Ri in IVCs. We also see 
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FIG. 4: Temperature evolution of Ic^ Iri and Ir2 for (a) an 
unshunted weak link (WL) and (b) the same WL with a par¬ 
allel shunt resistor, showing the reduction of crossover tem¬ 
perature in shunted case. The symbols are the data points. 
The continuous lines are the fits given by (in mA and K) 
Icin) = 0.18(7.2 - n) and Ir 2 (Tb) = 0.32(Tc2 - for 

both the devices. 


large fluctuations in voltage close to and for currents 
near in both the devices, see Fig. ie) and|^d). From 
the time-series data we find a bistable telegraphic-like 
voltage signal in this regime. Thus time averaged volt¬ 
ages in the IVCs show significant fluctuations. Close to 
the boundary of the bistable regime, more sensitivity to 
noise is indeed expected. 

We have also studied a shunted /i-SQUID device, with 
the same shunt geometry and resistance Rs value as in 
the single WLs discussed above. Although the critical 
current Ic is smaller in this SQUID due to the reduced 
width (< 50nm) of the weak links. The SQUID loop pat¬ 
tern is otherwise identical to our earlier work [8]. Figure 

shows the IVCs of the shunted /i-SQUID at different 
temperatures, which is found to be non-hysteretic down 
to 1.3 K temperature, establishing the role of the shunt in 
widening the non-hysteretic temperature range for both 
WLs and /i-SQUIDs. SQUID oscillations are clearly ob¬ 
served, see Fig. inset. 

Finally, let us discuss how we could further increase 
the re-trapping current and hence expand the hysteresis- 
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FIG. 5: IVGs of a shunted /i-SQUID at 1.3, 4.5, and 6.7 K 
showing non-hysteretic characteristics. The inset shows the 
SQUID oscillations at similar temperatures biased close to the 
critical current namely 0.37 mA for 1.3 K, 0.22 mA for 4.5K 
and 0.05 mA for 6.7 K. The vertical scale at 6.7 K has been 
enhanced by a factor of 5 for better clarity in large scales. 


free temperature range. Using a lower Rs value will in¬ 
crease Iri further and widen the hysteresis-free temper¬ 
ature range for a given WL device. Nevertheless, this 
will also reduce the overall normal resistance and result 
in a lower voltage signal to be measured. The same can 
also be achieved by using a smaller Ic WL with same Rg 
value. We have verified this claim in another shunted 
WL device, with same Rs value and smaller U, showing 
a Ths below 4.2 K. In any case, the shunt-resistor has to 
be kept close enough to the WL, so as to avoid relaxation 
oscillations, but not too close to cause heat or electron 
sharing between WL and the shunt, which can affect the 
WL superconductivity. 


CONCLUSIONS 

In summary, we have demonstrated a significant im¬ 
provement of the hysteretic behavior of superconducting 
WLs and /i-SQUIDs, using a parallel resistive shunt in 
close vicinity to the WL. As a result of the shunt, the 
hysteresis-free temperature range is wider. Our results 
can help to further develop WL-based non-hysteretic de¬ 
vices such as SQUIDs. 
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FIG. 6: (a) Schematic of the device showing the hot-spot in 
narrow leads, while the short weak link of length I act as a 
power source, (b) As/Ao variation with Rs as found from 
the generalized thermal model (symbols) and simple thermal 
model (thin continuous lines). These plots are shown for two 
different lengths, L = 3lt (red curve) and L = 9lt (blue curve), 
of the narrow leads. The inset in (b) shows the normalized 
IVGs for L = 3lt with Rs values of oc, 100 (marked with the 
arrow), 50, 35.6, 20 and 10 fl. Here, v{xo) = V{xo)/2IoRi. 
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SUPPLEMENTARY MATERIAL 

HOT SPOT MODEL FOR THE SHUNTED WEAK 
LINKS 

We discuss here a more elaborate quasi-static thermal 
model for our devices. This is a special case of the model 
discussed by Skoepol et. al. m- In our devices, the 
WL is quite short, as compared to thermal length, and 
the normal metal-superconductor (NS) interface occurs 
in the adjacent narrow leads which loose heat to the bath 
via the interface with the substrate. The heat conduc¬ 
tion, together with the interface heat loss, occurs in the 
narrow leads (of length A, width w and thickness t) that 
join the wide leads as shown in Fig. [^a). We assume that 
the temperature where the narrow lead meets the wide 
one is since the wide-lead is much wider than the nar- 
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row one. We treat the WL as a simple source of PRwl 
heat when it is resistive. The unsustainability of NS in¬ 
terface in the narrow lead determines the magnitude of 
the re-trapping current. Near the NS interface, the tem¬ 
perature is close to Tc and thus the thermal conductivity 
(K) does not have much variation near this interface. K 
is then considered as a constant. Solving the heat flow 
equations for the narrow leads, we get the relation be¬ 
tween bias current and the NS interface location, xq, as 
I = lQi{xo) with 


Io = 


lawH{Tc - Ti,) 
Pn 


( 2 ) 


and 


i{x) = 


sinh 1 

) [1 + coth 1 

(4 

j coth 1 

{it It) 

q + sinh ( 


1 


(3) 


Here, pn is the the normal state resistivity of the nar¬ 
row lead. It = ^jKt/a is the thermal healing length, 
which is defined as the length scale over which the heat 
(Qo = walt{Tc — Ti))) is transferred to the substrate. 
Here the coefficient a is characteristic of the interface. 
q represents the ratio of Rwl to the resistance of 2 lt 
length of the narrow lead. The voltage across the WL 
for a given xq is: V{xo) = Ioi{xo)[4:Ri{xo/L) + Rwl]- 
In the presence of Rs the bias current is given by / = 
Ioi{xo) + [V{xo)/Rs], which gives a different xq value, 
while the voltage remains the same. The inset of Fig. 
ib) shows the calculated IVCs for different Rs values 
with a current minimum at irs = Irs/^o^ which de¬ 
fines the normalized re-trapping current. Further, irs 
increases with reducing Rg value, see Fig. [^b). We find 
that irs value mainly depends on the narrow-lead and 
shunt resistance values but it does not depend much on 
Rwl- Iq follows a square root dependence on (Tc — T^) 
implying the same for Irs- The increment in re-trapping 
current by the parallel shunt resistance is in agreement 
with our model with Rg = 35.6 Rn = 2Ri = 56 Q and 
Rwl = 9 

Fig. [^b) shows the enhancement of retrapping current 
for two different lengths of the narrow lead. For compar¬ 
ison, we have also plotted results from the simple model 


(Eq. 1 of the main manuscript). Since the length of nar¬ 
row leads (L=3 It) in our devices is comparable to It^ we 
are still in a regime, where conduction dominates over 
surface heat loss. Thus both the models give similar en¬ 
hancement in retrapping current. With longer length of 
the narrow leads, the simple model deviates significantly 
from the detailed model as seen in Fig. [^b) for T = 9/f 



FIG. 7: IVC at higher bias currents (as an extension to Fig. 3 
of the main text) at 1.3 K temperature, for shunted WL. This 
shows the different transitions corresponding to the thermal 
runaway in various portions of the device. 


IVC OF THE SHUNTED DEVICE AT HIGH BIAS 
CURRENTS 


Fig. 0 shows IVCs at higher bias currents, where we 
see several jumps in voltage occurring at different cur¬ 
rents. At 1.3 K these currents are: Iri = 0.45 mA, 1^2 = 
0.78 mA, Ic = 0.85 mA, Rs = 1.80 mA, Ra = 2.08 mA, 
and Rs = 2.18 mA. We see several kinks between 1.82 
mA and 1.43 mA while ramping down the current. The 
jumps at Jc, /r 3 , R 4 and Rs are seen, when we ramp the 
current up, while the others occur during current ramp 
down. These transitions have been seen reproducibly in 
three different devices and are attributed to the thermal 
instability of N-S interface in different superconducting 
portions. The magnitude of R 2 agrees with the insta¬ 
bility current in 2 pm wide leads [8]. All the five re¬ 
trapping currents follow similar square-root dependence 
on (Tc — R) as expected for the instability currents [8]. 









